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Abstract 

To take into account the temperature fluctuation in the Planck distribution, we calculate con- 
volution integral with several probability distributions. Using these formula as well the Planck 
distribution and a formula in the non-extensive statistics, we analyze the data measured by the 
Cosmic Background Explorer (COBE). Our analysis reveals that the derivation from the Planck 
distribution is estimated as \q— 1| = 4.4 x 10 -5 , where q means the magnitude of the non-extensivity 
or the temperature fluctuation, provided that the dimensionless chemical potential proposed by 
Zeldovich and Sunyaev exists. Comparisons of new formulas and the Planck distribution including 
the Sunyaev-Zeldovich (S-Z) effect are made. 

PACS numbers: 05.10.Gg, 05.30. Jp, 98.70.Vc 
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I. INTRODUCTION 



One of interesting subjects in thermodynamics is relating to the non-extensive statistics 



Several years ago, authors of Refs. 



.1. 



showed that formulas in the non-extensive 



statistics, i.e., the temperature fluctuations, are calculated by the convolution integral with 
the gamma distribution 4|] and a calculation with a distribution described by exp(— 2|w| a ), 



where a is fractional number 
the Planck distribution 



.5] 



In this report, we calculate the temperature fluctuation in 



D 



Planck 



C R v 3 



(1) 



where Cb = 87rh/c 3 , u — hp and /3 = 1/ksT. \x denotes the dimensionless chemical potential 
introduced by Zeldovich and Sunyaev p. Equation (JTJ) is applied to the data measured by 







the Cosmic Background Explorer (COBE) 

The formula in the non-extensive statistics is obtained in Ref. 



12] 



as 



x \ 1 + (1 — q)x 



-- Dpi anck ((3, v)\\-e 
1 + e~ x x 1 + 3e 



-asi(g-l) 



(2) 



.l-e~ x 2(l-e" x ) 

where x = j3u. (It should be noticed that calculation in Ref. |2| is done without \i. The 



value of x 2 is 108 and the degree of freedom is 31 

In the next section, we calculate the convolution integrals using the several probability 



distributions l| 10, In the 3rd section we analyze the data reported in Ref. [12J, [l3j in 
terms of Eqs. (0) and (J2J) with /i, as well as new formulas. In the final section, concluding 
remarks are presented. We look for the origin of the temperature fluctuation. 



II. CONVOLUTION INTEGRALS OF EQ. © WITH PROBABILITY DISTRI- 
BUTIONS 



First of all, according to Refs. 
ature fluctuation as 



we consider the gamma distribution for the temper- 

P{<x,PM = f^t(^) >-V/^ (3) 

The formula, Eqs. Q, is also named the x 2 distribution. The convolution integral is given as 
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C jB z/ 3 e~ nA1 ~ A)a;n+ 5' 9 o w2?l2(9 ~ 1)+0((9 ~ 1)2) , (4) 



n=l 



where a = and (g-l)<l are assumed. 

The following generalized x 2 distribution is known as the non-central yf. distribution or 

1QQ. 



the generalized Glauker-Lachs formula and Perina-McGill formula 

a(l — p) a (3 
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where p = (nth)/ (ntotai) and J denotes the modified Bessel function, (nth) does the thermal 
boson. Eq. (jJJJ) is the solution of the following Fokken-Planck equation 



d* 7 



(6) 



where a and 77 are parameters. This is named the stochastic process for birth and death 
process with the immigration. The convolution integral of Eqs. (JTJ) and © is expressed as 



1 
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where 



M-) 0--p) 



a 

P(3o 



+ um 



As p -> 1 and a = ^ > 1, Eq. © reduces to Eq. (gj). 
As the third case, we consider the Gaussian distribution. 

1 r (P~Po) 2 



P(V, P, Po) 
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where V 2 denotes the variance. This distribution is relating to the solution of the Ornstein 
Uhlenbeck stochastic process, 



dP d „„ a 2 d 2 P 

= r Qp J 

dt dp ^ 2 dp 2 



(9) 
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FIG. 1: Analysis of data by the COBE in terms of Eqs. |JT|). The data sets are obtained from 



T. C. Mather and D. J. Fixsen. See Refs 



13j ]. We use a method by D. Seibert 



141 ] in calculation 



of error bar. The diagonal components are shown. 



From the convolution integral of Eqs. (0) and (jHJ, we obtain the following formula 



F (Gmss) (A, v, a, ij,V) = J2 C B u 3 e- nfl e 



(10) 



n=l 



When V 2 is expressed by /3g(g — 1), i.e., the temperature fluctuation, Eq. (fTUj) reduces to 
Eq. © 

III. ANALYSIS OF DATA BY COBE BY MEANS OF EQS. ©, ©, ©, © AND 



Using the CERN-MINUIT program, we can estimate values of parameters contained 
in formulas give n in previous sections. In Fig. Q we have confirmed the result reported 
in Refs. [12I Our results are shown in Table HI As seen in Table HI values of x 2 ' s 
for Eqs. (J2J and (@J) become to be somewhat better than that of the Planck distribution 
(Eq. ((H)), because of the parameter q. The magnitude of (q — 1) denotes the degree of the 
deviation from the Planck distribution including /1, i.e. Eq. ( |Tj) . 
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TABLE I: Analysis of data by the COBE by means of Eqs. ©,©,©, © and flU). 



formula Tq = 1//3q 1 X 2 / n -d-f M P 



Eq. © 


2.726 ± 0.001 




141/41 


(fixed) 




(Planck) 


2.72b ± 0.001 




O O / A f\ 

82/40 


(—3.32 ± 0.43) x 10 


-4 


Eq. ® 


2.726 ± 0.001 


(2.00 ± 0.07) x 10~ 5 


167/40 


(fixed) 






2.725 ± 0.001 


(4.57 ± 1.60) x 10 


t A /on 

74/39 


(—5.15 ± 0.77) x 10 


-4 


Eq (131) 


2 726 ± 001 


(2 00 ± 06) x 10~ 5 


169/40 


(fixed) 




(gamma) 


2.725 ±0.001 


(4.40 ± 1.50) x 10" 5 


74/39 


(-5.16 ±0.76) x 10" 


-4 


Eq. © 


2.726 ± 0.001 


(-3.99 ±0.86) x 10~ 5 


119/39 


(fixed) 


0.9995 ± 0.6842 


(NCC) 


2.725 ±0.001 


(4.40 ± 1.62) x 10~ 5 


74/38 


(-5.16 ±0.74) x 10" 


4 0.9995 ± 0.6641 


Eq. (ID) 


2.726 ± 0.001 


(2.00 ± 0.06) x 10~ 5 


169/40 


(fixed) 




(Gauss) 


2.725 ±0.001 


(4.40 ± 1.47) x 10~ 5 


74/39 


(-5.16 ±0.74) x 10- 


-4 



IV. CONCLUDING REMARKS 

We have calculated Eqs. (j3J), (JTj) and (fTUj) and applied them to the data measured by 



the COBE 



12, 



3 



. The following remarks are obtained from our analysis. 



1) It is shown that the dimensionless chemical potential /i is playing an important role 
in Eq. (JTJ), because of values of x 2 's- 

2) The deviations from the Planck distribution including the chemical potential for the 
cosmic microwave background (CMB) photon is estimated by \q — 1| = 4.4 x 10 -5 . 

Both approaches, the non-extensive statistics and the convolution integrals, show the 
same magnitude. 

3) As (q — 1) <C 1 and V 2 = Po(q — l) are assumed, it is difficult to distinguish the gamma 
distribution and the Gaussian distribution for the temperature fluctuation. 

4) From p w 1 (in spite of large error bars) in Eq. ((7j) in Table |H it is known that the 
CMB photon is the pure thermal photon. 

5) To look for the origin of the magnitude of the non-extensivity or the temperature fluc- 
tuation, we analyze the data by the following formula including the Sunyaev-Zeldovich 



TABLE II: Analysis of data by the COBE in terms of Eq. (JTTJ) . 



formula 


To = 1/A) 


y 


X 2 /nA.f 




Eq. CEU 


2.726 ±0.001 


(-1.99 ±0.43) x 10~ 5 


119/41 


(fixed) 




2.725 ±0.001 


(2.20 ±0.75) x 10~ 5 


74/40 


(-5.16 ±0.76) x 10~ 4 



(S-Z) effect IslQ, 



F (S-Z effect)^ ^ ^ y) = Dpi anck (f3, V, /i) + C B U 3 



xye 



i x 
x coth — 



where x = (3 u + \i. y is named the y parameter in the Compton scattering, 

y=[ dln e a T ^\ , (12) 
J m e c 2 

where I, n e , ot and T e are the size of region with high temperature in the cosmos, 
the number density of electrons, the cross section of Thomson scattering and the 
temperature of electron, respectively. Our estimated values are shown in Table |H] 
The magnitude of y parameter y « (2.2 ± 0.75) x 10~ 5 is almost the same as that of 
\q — 1| in Table E This coincidence suggests that Eqs. (jlj), and (flT!|) work well. 

We summarize the following concerning the item 5): The deviations from the Planck 
distribution for the CMB photon are seen through the non-extensivity or the temperature 
fluctuation. Moreover, the physical origin is attributed to the Compton scattering between 
electrons in the galaxies and CMB photons 0, Q|. The behavior of the second term of 
Eq. (fTT|) is shown in Fig.|21 The crossing point on the horizontal line (frequency) is determined 
by x coth f = 4. Differences between Eqs. © and and the Planck distribution with fi are 
also shown in Fig. |2l The small quantity, (0.1-0.05)/400, is reflecting the physical meaning 
of the inverse Compton scattering in the galaxies. It is observed through Eqs. and 
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FIG. 2: Correction term to the Planck distribution. Planck distributions are subtracted at data 
points, (a) The second term of Eq. 1)11(1 . i.e., Sunyaev-Zeldovich (S-Z) effect, (b) Difference 
between Eq. © and the Planck distribution with /i. (c) The same as (b) but Eq. @ and the 
Planck distribution with /i. 
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